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ABSTRACT 

The recent discovery by M. Konacki of a "hot Jupiter" in the hierarchical triple star system HD 188753 
challenges established theories of giant-planet formation. If the orbital geometry of the triple has not changed 
since the birth of the planet, then a disk around the planetary host star would probably have been too compact 
and too hot for a Jovian planet to form by the core-accretion model or gravitational collapse. This paradox is 
resolved if the star was initially either single or had a much more distant companion. It is suggested here that a 
close multi-star dynamical encounter transformed this initial state into the observed triple, an idea that follows 
naturally if HD 188753 formed in a moderately dense stellar system — perhaps an open cluster — that has since 
dissolved. Three distinct types of encounters are investigated. The most robust scenario involves an initially 
single planetary host star that changes places with the outlying member of a pre-existing hierarchical triple. 
Subject headings: open clusters and associations: general — planetary systems — stars: individual (HD 
188753) — stellar dynamics 



1. INTRODUCTION 

iKonackil J2005I) recently reported the detection of a giant 
planet in star system HD 188753. The planet orbits a Sun-like 
star every 3.35 days and has roughly the mass of Jupiter. Al- 
though this orbit is extremely tight, it is unremarkable among 
~30 known "hot Jupiters" with periods of 2.5-10 days around 
solitary stars 1 . What distinguishes the new discovery is that 
the planetary host star (A) is bound to a companion (B) that 
is itself a compact binary; HD 188753 is thus a hierarchical 
triple star. Star A has a mass of M A = 1.06 ± 0.07 M Q , and 
binary star B has constituent masses of Mbi = 0.96 ± 0.05 M Q 
and Mb2 = 0.67 ± 0.05 M Q . Orbit AB has a semimajor axis of 
aab = 12.3 ±0.04 AU and an eccentricity of e AB = 0.50±0.05, 
and B has parameters «b = 0.67 ± 0.01 AU and e% = 0.10 ± 
0.03. While the notion of a planet in a triple star is quite 
provocative, the most fascinating aspect of HD 188753 is the 
close orbit of AB, which seems to challenge existing models 
of giant-planet formation. 

Planets form in a disk of gas and dust surrounding a nascent 
star. Two theories exist for the genesis of Jovian planets (see 
Bodenheimer & Lin 2002 and references therein). In the core- 
accretion paradigm, a solid core of ~10 Earth masses forms at 
disk radii of >2-3 AU, where the temperature is low enough 
for ices to condense. The core then accretes a gaseous enve- 
lope. The alternative picture involves the gravitational frag- 
mentation of the disk into gaseous protoplanets at radii of ~5- 
10 AU. When addressing hot Jupiters, each theory requires 
some mechanis m to drive the migration of the plan et to small 
radii (e.g., Lin, Bodenheimer, & Richardson 1996). 

Complexities and uncertainties inherent in these two mod- 
els are exacerbated when the young star is in a binary. The 
perturbing gravitational field of the binary c ompani on acts 
to (1) tidally truncate the disk (e.g., [A rtvmowicz & Lubow 
Il994t iPichardo. Soarke. & Aguilarl 1 20051) and (2) stir and 
heat the disk (e.g.. iGoodmanll 19931 iNelsonl 120001) . A disk 
truncation radius of <3 AU may p reclude the embryonic stage 
in the core-accretion scenario (lJang-Condellll200'5l) . Even if 
the truncation radius is ~10 AU, the disk may be too hot to 
allow either core formation or unstable gravitational collapse 

^http : / /vo . obspm . f r/exoplanetes/encyclo/encycl ■ html| 



JNelsonl2000l) . 

The source of the conundrum in HD 188753 is that the or- 
bital parameters of AB imply that a disk around star A would 
be truncated at a radius of only ~1 AU. It seems very unlikely 
that a giant planet could have grown in such an environment. 
However, this dilemma can be avoided if we abandon the tacit 
assumption that the system has maintained its current config- 
uration since the time when the stars and planet formed. 

Suppose that star A and the planet have always been paired, 
and that the planet was allowed to form by one of the two 
models sketched above and then migrate to small radii. It fol- 
lows that, at the time of planet formation, star A was either 
single or in a much wider binary than at present. The simplest 
way to transform either of these states into HD 188753 is via 
a strong dynamical interaction, such as a scattering encounter 
between a single and a triple star, or two binaries. This could 
only have occurred in a moderately dense stellar environment. 
Different modes of multiple-star formation and probabilistic 
arguments regarding the birthplace of HD 188753 are dis- 
cussed in § Three plausible dynamical histories for HD 
188753 are introduced in §|3] In §0] numerical scattering ex- 
periments are used to gauge the likelihood of each scenario. 
The broader relevance of this work is given in §|5] 

2. CLUSTER PROPERTIES 

Stars rarely form in isolation, and so it is quite likely that 
HD 188753 was born in some flavor of stellar system, either 
a loose, transient group of stars or a bona fide open cluster 2 . 
The salient features of groups and open clusters are listed be- 
low (see Adams & Myers 2001 for a thorough discussion). 
Once the relevant kinematical parameters are in hand, it is a 
simple matter to assess the most probable environmental con- 
ditions for a dynamical origin of HD 188753. 

Most stars are born in small groups with N ~ 10— 10 2 mem- 
bers, stellar densities of n < 100pc" 3 , velocity dispersions 
of a < lkms -1 , and evaporation time scales of T < lOMyr. 
There is insufficient data on these short-lived groups to de- 
velop a reliable statistical picture of their properties and stel- 
lar content. Roughly 10% of all stars form in larger sys- 

2 The triple did not form in an old Galactic globula r cluster, becau se, 
among other reasons, star A is too young and massive (see Konacki 2005). 




tems (N ~ 10 2 -10 4 ) that would rightfully be called open clus- 
ters, for which there is a wealth of observational data. Open 
clusters tend to have higher densities and somewhat higher 
velocity dispersions than groups, and have longer lifetimes. 
While some clusters are as old as ~10 10 yr, the median age 
is ^10 8 yr. All but a few percent of open clusters dissolve in 
< 10 9 yr, due to the dynamical effects of internal evolution, the 
Galactic tidal field, and close passages by molecular clouds 
(e.g.. lWielenll985l) . In open clusters that have been surveyed 
for stellar multiplicity, the estimated binary fraction is ^50% 
(e.g.,|Bica & Bonatto 200 3) ■ Although few triples are k nown 
in clusters (e.g.,[Mermilliod, Duquennov, & Mavo ilT994l) . the 
expected fraction is ~5 -10%, based on the s tatistics of field 
stars near the Sun (e.g., Dua uennov & Mavorl 19911) . 

A dynamical interaction capable of producing HD 188753 
involves at least one binary (or hierarchical triple) with an 
orbital separation of a ~ 10 AU (see §[3j. A binary of this 
size can be strongly perturbed by a third star if the distance 
of closest approach between the star and binary center of 
mass is <a. For stellar speeds of a ~ lkms" 1 , we see that 
GMq I a 1 ~ 10 3 AU ^> a. Therefore, the close-encounter cross 
section is dominated by gravitational focusing and takes the 
characteristic value S ~ 2naGM j a 2 , where M is the total 
mass of the few-body system. The probability that the binary 
has such an interaction in the lifetime of the stellar system is 



P ~ n£oT : 
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where n 2 = n/100pc" 3 , o 10 = a/lOAU, M 3 = M/3M©, a\ = 
er/lkms" 1 , and Tj = r/10 7 yr. Suppose that the number of 
similar binaries in any stellar system is a fixed fraction of the 
total number of stars. The net probability for a strong few- 
body interaction to occur in a given system would then scale 
as NP. A t ypical open cluster has N ~ 10 3 and Tj ~ 10 (e.g., 
lFriel l995). making the characteristic value of NP for a clus- 
ter ^100 times larger than for a smaller group. Even when 
we multiply by the additional factor of 0.1, the fraction of 
stars formed in open clusters, it is still very probable that HD 
188753 originated in a system with ^10 3 stars that lived for 
>10 8 yr and has long since evaporated. 

3. DYNAMICAL HISTORY OF HD 188753 

It may seem a daunting task to investigate all possible ways 
of forming HD 188753 in a cluster rich in multiple-star sys- 
tems. However, three distinct scenarios stand out as being 



the most likely (see Fig^for graphical depictions): (I) hi- 
erarchy AB existed at the time of planet formation, but with 
a much larger periastron separation, which was subsequently 
reduced to its present value in a strong encounter with a single 
star; (II) star A was born with a different, single companion 
in a wide orbit that was later exchanged for binary star B; 
(III) star A was initially single before being exchanged into 
a pre-existing hierarchical triple containing B. Each channel 
involves a catalyst star (C) that ultimately escapes. In all of 
these interactions, the finite sizes of the planetary orbit, binary 
B, and star C are neglected, an idea that is exploited in §0] 

In scenarios I and II, star A is in a binary with B or C at 
the time of planet formation. The semimajor axis a and ec- 
centricity e of this binary are constrained by the demand that 
the planet formed by the core-accretion process or gravita- 
tional collapse at a radius of >3 AU in a disk around star A. 
iPichardo. Sparke. & A guilarl 112005) estimate the disk trunca- 
tion radius to be R, = \a{ \ -e) 12 , where A, a function of the 
stellar masses, decreases from ~l/3 to ~l/5 as the com- 
panion mass increases from O.2M to 2M Q . It is not clear 
how large the disk and binary orbit must be to allow giant- 
planet formation to proceed much as it would around a single 
star. The simple constraint adopted in the numerical study 
of the next section is that R, > 10AU, which implies that 
a(l-e) 12 > 40 AU in scenarios I and II. 

In scenario III, single star A encounters triple BC, where B, 
a tight binary with separation a B , orbits C with a semimajor 
axis of abc and an eccentricity of <?bc- It seems likely that BC 
will have long-term stability against disintegration, for which 
an approximate condition is (Mardling & Aarseth 2001) 



flBc(l-<?Bc) > 2.8a B 



M B c (1+cbc) 



Mb VI ~<?bc 



2/5 



(2) 



where M BC = M B +M C . When a B = 0.67 AU, M B = 1 .63 M , 
and Mc = 1M Q , the stability condition is abc(1 _ <?bc) 1,2 > 
2.3-3 AU for <?bc = 0-1. Note that scenario III is much less 
constrained than I and II. It is shown in the next section that 
scenario III is the favored channel for producing HD 188753. 

4. NUMERICAL INVESTIGATION 

Simulations have been conducted to demonstrate how ef- 
fective each of the three scenarios is in producing HD 188753. 
Instead of following all 5 objects — 4 stars and 1 planet — the 
problem has been reduced to scattering encounters between 
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a AB (AU) 

FIG. 2. — Each panel shows the semimajor axis and eccentricity of every 
pair AB that is produced in 10 5 scattering experiments; the number of dots is 
given in parentheses. Red stars mark the orbital parameters of HD 188753. 
The curves are the stability boundaries for the triple (see eq. |2*1). 



binaries and single stars. This idealization is justified when 
no pair in the set {A,B,C} comes sufficiently close to induce 
strong tidal perturbations. The adopted radii of the objects 
are the orbital separation of the planet from A (0.045 AU), 
the semimajor axis of B (0.67 AU), and the radius of a main- 
sequence star for C, approximated as 0.005(Mc/M Q ) AU. A 
scattering calculation is terminated if the separation of any 
pair is less than three times the sum of their radii; the factor 
of three is a somewhat arbitrary, but conservative, choice. In 
practice, a minor fraction of the experiments end this way. 

Initial conditions for the scattering encounters are gener- 
ated by Monte Carlo methods from astrophysically motivated 
probability distributions. The orientatio ns and orbital phases 
of the incident binaries are chosen as in Hut & Bahcall ( 1983). 
Semimajor axes are drawn from a logarithmically flat distri- 
bution, p(a) oc a~ l (a = 3-100 AU), ins pired by binary statis- 
tics i n the Solar neighborhood (e.g., Duauennov & Mavor 
119911) . A linear distribution, p{e) oc e (e = 0-1), is adopted 
for the initial eccentricities. Impact parameters are distributed 
uniformly in b 2 , with a maximum value large enough to 
allow all outcomes of interest. Relative speeds at infin- 
ity, v, are picked from a Maxwellian distribution, p(v) oc 
v 2 exp(-v 2 /2<7 2 ). There is no distribution of catalyst masses, 
Mc, that is suitable for all three scenarios. As a concrete, but 
rather arbitrary, choice, Mq is drawn from p(Mc) oc M^ 1 over 
the range of 0.2-2M Q . 



Once the parameters of the incident binary have been cho- 
sen, they are checked for consistency with a given scenario. In 
scenarios I and II, R, > 10AU must be satisfied (by assump- 
tion) for the incident binary containing A. Scenario III only 
requires that the incident triple is stable (eq. |3). Each sys- 
tem that passes this initial test is followed through the scat- 
tering interaction. The code used to automatically generate 
and termin ate each encounter is similar in logic to that de- 
scribed in Hut & Bahcall (1983). The equations of motion 
ar e integrated with the time-tr ans forme d leapfrog algorithm 
of Mikkola & Tanikawal(ll999l) and lPreto & TremaineM 19991) . 

Every time a bound pair AB is left as the final product of an 
experiment, the orbital parameters are recorded and the sys- 
tem is tested for triple stability; unstable systems are cut. Fig- 
ure |3 shows the final results of these simulations. A total of 
10 5 scattering experiments were computed for each scenario. 
It is immediately apparent in Fig. |2] that scenarios I and II 
have considerable difficulty in producing systems that resem- 
ble HD 188753, while the corresponding efficiency for sce- 
nario III is quite high. If the minimum disk truncation radius 
were reduced to, e.g., 5 AU, scenarios I and II would be some- 
what more profitable. Nonetheless, scenario III is clearly the 
most robust of the three channels, simply because there are no 
special hurdles for Jovian-planet formation. 

5. DISCUSSION 

The obstacles to forming a giant planet in HD 188753 can 
be overcome if the triple originated in a moderately dense stel- 
lar system and has a nontrivial dynamical history. A simi- 
lar picture might apply to other close binaries or triples that 
host a giant planet. Three systems are noteworthy in this re- 
gard, each with an orbital separation of ~20AU: 7 Cephei 
jHatzes et aTll2003l): HP 41004 (a triple containing a brown 
dwarf: IZu7ker*etal.l2004): Gliese 86 ( has a white-dwarf sec- 
ondary; Mugrauer& Neuhauser 2005). All of these systems 
could threaten conventional theories of Jovian-planet forma- 
tion, but the compact orbits may have resulted from exchange 
encounters or other types of strong dynamical interactions in 
a stellar group or open cluster. 

Stellar dynamics can play critical roles in the planet- 
formation process and in generating variety among single, 
binary, and triple stars that host planets. More than half 
of all field stars are in binaries and higher-order multiples, 
and most of these stars are born in systems with >10 
members. Thus, it can be argued that stellar dynamical 
considerations are fundamental to a complete understand- 
ing of how planets are created and evolve dynamically, 
and in deciding how to search for these worlds. Earlier 
studies of plan et formation and dynamics in star clus- 
ters (e.g.. |de la Fu ente M arcos & de la Fuente Marcos 
1 19971 iLaughlin & Adamsl 119981 iBonnell et alJ 120011 
Hurlev & Shara 120021) highlight the complexities and 
richness of these problems. However, much work remains to 
be done. Investigations based on scattering experiments, like 
those conducted here, should be expanded to include binary- 
binary, single-triple, and binary-triple stellar interactions. 
Alternatively, detailed A^-body simulations should be used 
to study the fates of planets in systems with ^10-10 4 stars, 
covering a wide range of initial conditions and assumptions 
regarding the statistics of binary and triple stars. 

A recent, inde penden t study by 

Portegies Zwart & McMillan (2005) uses similar meth- 
ods of calculation to those presented here and comes to 
essentially the same conclusions regarding the history of HD 
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188753. 
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